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ABSTRACT 

The exact ionization state of Be^ in the solar core is crucial for the precise pre- 
diction of the solar neutrino flux. We therefore examine the effect of pressure 
ionization on the ionization state of Be*" and all elements with 12 > Z > 4. We show 
that under the conditions prevailing in the solar core, one has to consider the effect 
of the nearest neighbor on the electronic structure of a given ion. To this goal, we 
first solve the Schrodinger and then the Kohn-Sham equations for an ion immersed in 
a dense plasma under conditions for which the mean interparticle distance is smaller 
than the Debye radius. The question of which boundary conditions should be imposed 
on the wave function is discussed, examined and found to be crucial. 

Contrary to previous estimates showing that Beryllium is partially ionized, we 
find that it is fully ionized. As a consequence, the predicted rate of the Be^ + e" 
reaction is reduced by 20-30%, depending on the exact solar model. Since Be*" is a 
trace element, its total production is controlled by the unchanged He^ -I- He^ reaction 
rate, and its destruction is determined by the rate of electron capture. As the latter rate 
decreases when the Beryllium is fully ionized (relative to the case of partially ionized 
Be) , the estimate for the abundance of Be^ increases and with it the B* neutrino flux. 
The increase in <^y(B*) is by about 20 — 30%. The neutrino flux due to Be^ electron 
capture remains effectively unchanged because the change in the rate is compensated 
for by a change in the abundance. Hence the prediction for the ratio of ^j/(B®)/(/)^(Be^) 
changes as well. 

Key words: Sun:interior, equation of state, atomic processes, plasmas 



1 INTRODUCTION 



Classical calculations of solar model s assume that all species 



are fu lly ionized above ~ 10 K (e.g. Bahcall & Pinsonneault 
1992). The main reason probably being the saving of com- 



puter time, because detailed ionization calculations are very 
CPU demanding. Alternatively, one can use tables for the 
equations of state which are calculated to high accuracy. 
However, one then has to interpolate for the relative abun- 
dances which change continuously once diffusion takes place. 
Indeed, from the point of view of the total gas pressure and 
other thermodynamic quantities, the partial (or complete) 
ionization of heavy species like C,N,0,Mg or even Fe, af- 
fect the number of free electrons at temperatures above a few 
million degrees at a relative level of about 10"'^, depending 
on the exact mass fraction of the heavy elements. Conse- 
quently, the total pressure and speed of sound are affected 
at the same relative level of accuracy. 



Iben Kalata & Schwartz (1967, hereafter IKS67) exam- 
ined the ionization state of Be in the solar core, and con- 
cluded that its K-shell electrons are partially bound (with 
a population level of about 30% depending on the exact 
location in the core). This fact significantly affects the pre- 
dicted neutrino flux from the Sun. The most important 
channel for the destruction of Be^ in the Sun is via electron 
capture, of which most are free electrons. However, if the 
Be^ ion has some bound electrons then the rate of electron 
capture is enhanced, and with it the generated Be^ elec- 
tron capture neutrino flux (by about 20-30% once averaged 
over the entire relevant region in the Sun). Thus, the exact 
occupation fraction of the Be^ K-shell is important for the 
accurate prediction of the solar neutrino flux, and the ratio 
<^^(B*)/<j!>^(Be^) in particular. To include the effects of the 
plasma, IKS67 assumed a Debye Hiickel (DH) potential and 
calculated the dependence of the ground state energy on the 
environmental conditions. 
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The problem of obtaining the ion ization state o f Beryl- 
hum in the Sun was later revisited by Johnson et al. (1992) 



The authors analyzed the validity of the Debye-Hiickel po- 
tential and found that the prerequisites for the validity are 
weakly violated. The authors claim that once the assump- 
tions for the validity of the DH are strongly violated, "ex- 
periments show that the DH fails dramatically". In partic- 
ular, we note the first point raised by the authors, namely 
the requirement to have many particles in a Debye sphere 
needed for the validity of the DH treatment. This require- 
ment implies that the interparticle distance is significantly 
smaller than the Debye radius. The authors solve for the 
Beryllium atom assuming a DH potential using three differ- 
ent methods (DH, self-consistent DH and Hartree) and find 
only small differences in t h e ioni zation compared to IKS67. 
Gruzinov & Bahcall (1997) discussed the ionization 



state of Beryllium in the Sun assuming mean field screening, 
the density matrix formulation, and a Monte Carlo method. 
However, all within the framework of a screened Coulomb 
potential. The authors also discussed the effects of fluctua- 
tions and found only minor effects. 

If indeed Beryllium is partially ionized in the solar core, 
namely, it keeps the K-shell electrons at least part of the 
time, then several additional consequences follow. These ef- 
fects were hitherto neglected in the prediction of the solar 
neutrino fluxes (predictions that assumed at the same time 
partial ionization of Be^). 

First, the screening of the nuclear reaction Be^ -I- p, 
which is the competing Be'^ destruction reaction, should be 
calculated using the proper effective charge of the Be ion. 
Currently it is assumed to always have a charge of -|-4e. This 
effect would decrease the B* neutrino flux. 

Second, a similar correction to the screening should ap- 
ply to the higher Z reactions of CNO-|-p, affecting in this way 
the (small) contribution of the CNO cycle to the total solar 
energy budget. This effect would suppress the CNO energy 
production because the effective charge of the ion would be 
smaller and hence the electrostatic screening energy would 
be smaller as well. 

Third, the exact point at which various ions become 
completely ionized affects the entropy density in the outer 
convective zone of the sun, and with it the solar structure. 

As the accurate prediction of the solar neutrino flux is so 
important, the purpose of this contribution is to re-examine 
the ionization state of the heavy species in the solar core, 
and in particular the ionization state of the trace element 
Be^. 

The question to which extent does the Be^ or any other 
heavy ion, retain its K-shell electrons is usually analyzed in 
two steps (IKS67). The first step is to apply the simple Saha 
equation assuming that the structure of the Be'^ atomic en- 
ergy levels is unaffected by the dense plasma. The second 
step is to account for the effect of the plasma on the en- 
ergy levels of the Be^ by assuming a smooth DH potential 
and calculating the energy levels under the DH potential. 
Once the new energy levels are known, the electron popula- 
tion in the levels can be re-evaluated using a Saha equation 
which incorporates the revised energy levels. This approach 
is justified only as long as the plasma effects are small per- 
turbations. 

As we shall show, the conditions in the core of the Sun 
are very peculiar and the number of particles inside a Debye 



sphere, Nd ~ few and not very large compared to unity. 
Hence, the necessary condition for the validity of the Debye 
Hiickel theory is not satisfied. Moreover, the conditions in 
the core of the Sun are such that the mean distance between 
ions in general, (r^) = (47rn/3)~^''^, and between a proton 
or an a-particle and a Beryllium ion in particular, is of the 
order of 2Rb{Z = 4) (the Bohr radius in a nucleus with 
charge Z — 4) and hence the picture of an ion with an elec- 
tronic shell inside a Debye-Hiickel potential, is not strictly 
valid. Here, n is the number density of ions while the index 
s in (rs) corresponds to a calculation employing spherical 
packing of ions. When the DH theory applies, it means that 
there are many ions inside a Debye radius and the mean dis- 
tance between ions is much smaller than the Debye radius. 
The electronic structure of the ions is then affected first and 
foremost by the close ion rather than by the Debye cloud 
and its large radius. This point, which is essential in this pa- 
per, will be discussed at length since this situation dictates 
a different boundary condition which subsequently leads to 
different energy levels and a different ionization state (under 
the same thermodynamic conditions). 

The paper is structured as follows. We first question and 
analyze the effective potential to be used under the condi- 
tions prevailing in the solar core. Then we repeat the two 
steps analysis of IKS67. We next proceed to examine the 
pressure ionization at T = assuming a Coulomb poten- 
tial. In view of the doubtful validity of the DH potential 
under the conditions prevailing in the core of the Sun, we 
repeat the calculation assuming the Schrodinger and later 
the Kohn-Sham equations. We find that Be^ is fully ionized 
at a lower density (and temperature) then previously calcu- 
lated. Finally, we examine the effect of the complete ioniza- 
tion of Be'^ on the predicted solar neutrino flux according to 
different sets of nuclear reaction cross sections. 



2 WHICH POTENTIAL? 

The classical calculation of the atomic energy levels and 
pressure ionization assumes a smooth (in time and space) 
potential. We flrst examine the assumption that any ap- 
plied potential can be assumed to be smooth for the specific 
problem of the structure of the electronic level of a given ion 
under the conditions prevailing in the core of the Sun. 

For the assumption that the potential is smooth to be 
valid, the fiuctuations of the plasma must be much faster 
than the motion of the electron in the bound orbit such that 
the average smooth value can be taken (for the calculation 
of the bound state) rather than the instantaneous one. We 
assume at the beginning that the plasma does not affect the 
energies of the bound levels. (The effect of the plasma is 
to make the energy levels shallower so the electrons would 
move even slower. Thus, the argument presented here is a 
conservative one.) 

The Bohr radius in a Hydrogen-like ion with a charge 
Z (in vacuum) is given by: 



u2 1 

"0 = = ^0.528 X 10"*cm. 

The classical velocity is: 

1/2 



2Ze- 



(1) 



(2) 
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and the period is 



P = 



(3) 



There are two sources for fluctuations in the potential, those 
caused by the protons and those caused by the electrons. The 
typical time scales of the fluctuations are 



{rs)m. 



1/2 



(4) 
(5) 



where (r^) = f /(47rnion/3)'^''^ is the mean interparticle dis- 
tance and Tiion is the number density of ions. Vth is the rele- 
vant thermal velocity. The number of particles in the Debye 
sphere is given by: 



Nn = ^Rl, 
where 



Rd = 



kT 



(6) 



(7) 



is the Debye radius for this mixture, rij is the number den- 
sity of specie j with charge Zj. The above expression for 
the Debye radius assumes that both the electrons and ions 
contribute to the supposedly DH potential. For simplicity 



we assume TV 



1/2 



3 and obtain that for n = 10 #/cm , 



T = 1.5 X lO^K and a pure Hydrogen plasma, 



P 



0.5 and 



0.01. 



(8) 



Note that when some of the ions are Helium ions (in the 
core of the present Sun about half the ions are He), Nd 
decreases even more. What is the implication of this result? 
The fluctuations due to the electrons are of about the same 
timescale as the period of the electrons in the K-shell in a 
single ion in vacuum, while the protons in the plasma have a 
much longer time scale. Hence, it is not justified to treat the 
contribution of the electrons to the DH potential (felt by the 
electron) as a smooth potential in time. On the other hand, 
since there are only few protons in the Debye radius their 
contribution to the potential is smooth in time but not in 
space and certainly not spherical. Additional argum ents that 
ques tion th e validity of the potential are given by Johnson 
et al, (H992| ). In what follows, we assume that all potentials 
are temporally smooth and spatially spherical. 



3 THE IONIZATION OF BERYLLIUM IN THE 
SUN 

3.1 The state of ionization ignoring plasma effects 

Next, we discuss the ionization state of Be^ in the core of 
the Sun - the classical way. If one adopts the Saha equa- 
tion, ignoring screening and the excited energy levels (thus 
including only the ground states in the partition functions), 
then the probabilities /i and /2 that one or two K-shell elec- 
trons are associated with any given Be^ nucleus are given by 



(IKS67) 

/i = 

/2 = 
where 



V 



l+rj + 0.25r?2 exp{-Sx/kT) ' 



0.2577 exp(-<5x/fcr)/i 



' 2-KmkT 



(9) 



(10) 



Here xi = 216.6 eV is the forth ionization potential of 
the Be^ atom, X2 ~ 153.1 eV is the third ionization potential 
of the Be^ atom, and Ax = Xi ~X2 ~ 63.5 eV. These values 
correspond to the limit of vanishing plasma density. Ue is 
the number density of the free electrons, most of which are 
contributed by Hydrogen and Helium and are independent 
of the state of trace elements like Beryllium. Thus, can 
be treated as fixed. 

The application of the Saha equation in the above form 
to the core of the Sun (p = 158 g cm"^, T = 1.57 x lO^K, 
X = 0.36 k Z = 0.02) yields fi = 0.320 and /2 = 0.038, 
implying that the Be^ keeps its last electron for about a 
third of the time. 

As the relevant ions are in a plasma, the traditional pro- 
cedure to correct for the plasma eff ect is to replace th e pure 
Coulomb potential with a DH one ( |^ogers et al. 197C ) . This 
is for example the procedure IKS67 evaluated the plasma 
corrections for the energy levels of Be^ in the solar core. 

The Saha equation in the above form ignores electron 
degeneracy, exchange effects and pressure ionization. The 
electron degeneracy introduces a small correction under the 
conditions prevailing in the Sun. As we will shortly demon- 
strate, exchange and pressure ionization are significantly 
more important. In what follows we do assume in spite of 
the previous reservations, a smooth static DH potential con- 
tributed by the electrons and ions. Moreover, we assume it 
to be relevant in a statistical sense only. 



3.2 Screened potential: taking into account the 
plasma effects 

After performing the above estimate, IKS67 turned to evalu- 
ate the ground state of the Z = 4 ion assuming a smooth DH 
screened potential in which both the protons and the elec- 
trons are taken into account. We ignored the questions raised 
in the previous section concerning the validity of the poten- 
tial for our particular purpose here (ionization in the core 
of the Sun), repeated their calculation and confirmed their 
results with r e spect to the Hydrogen like ion with Z — 4. 
Rogers et al. (197C) calculated the bound states of static 



screened coulomb potential and formulated their results in 
terms of the screening length. We also repeated their results 
for the Is state and the results are shown in fig. 0. Clearly, as 
the Debye length approaches the Bohr radius of ions with 
charge Z, there are no more bound states. The boundary 
condition on the wave function in this case is i/)(r — > cx) = 0. 

The calculation of the plasma effects on the triply ion- 
ized Be^ ion is more complicated because of the partial 
screening of the nucleus by the bound electron. To overcome 
this problem, we used the following approximate method. 
We looked for the eigenvalue in the low density limit and 
searched for the effective charge that will reproduce the mea- 
sured ionization potential of 153.1 eV. We found that this 
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Figure 1. The energy of the ground state for a screened Coulomb 
potential as a function of the screening length. Note the units of 
the energy (^^ERydberg) s-^id Debye length scale {Rd / Rb{Z)) 
where Rg (Z) is the Bohr radius of an ion with charge Z. 




Density(g cm"3) 

Figure 2. The run of the ground energy levels of a Be ion with 
one and two electrons as a function of density for a Debye Hiickel 
potential and ■!/>(r — > oo) = as the boundary condition. The tem- 
perature, which enters via the Debye radius, is taken as constant 
at 1.57 X lO^K. Note that the temperature in the Sun decreases 
with density and hence the temperature and with it the Debye ra- 
dius are overestimated here for densities lower than 150 g cm"-^. 
The correction for the accurate temperature is small. The pur- 
pose of the figure is to show the conditions in the solar core. The 
small arrow marks the density at the core. 




0.4 0.6 
m{r)/Msun 

Figure 3. The occupation numbers f\ and /2 of Be throughout 
the Sun as a function of the solar mass fraction. The broken lines 
are the results of IKS67. The continuous lines are the present 
results after incorporating the effective charge of the Be'^+ ion 
(Zoff = 3.354). The curves marked with fi{KS) and f2(KS) are 
the results assuming i/i' = on the cell boundary. These results 
are the actual run of the occupation numbers in the Sun. Also 
shown are the density and temperature (in units of the central 
temperature) in the present day Sun. The run of Ro/Rb and 
{rz)/RB are shown as well. T/Tc, p, Rd/Rb and {'^z)/Rb are 
all shown on the left axis. The occupation numbers arc shown on 
the right axis. 
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Figure 4. The run of the energy of the ground state as a func- 
tion of density assuming fixed temperature and composition. Also 
shown are i?£)/i?s(Z = 4) and {rs)/i?s(Z = 4). The arrow 
marks the density in the center of the Sun. 



charge is Z = 3.3544. We then repeated the calculation of 
the plasma effect on the bound states assuming this charge. 
The results are shown in fig. ^. We then used the new values 
for the ionization potential in the Saha equation to find the 
revised /i & /2. The comparison between the values used 
with and without the plasma correction are shown in fig. ^ 
where the actual run of the occupation numbers in the Sun 
is given. It is surprising that the differences in the ionization 
come out to be quite small. 



The particular results for the binding energy (calcu- 
lated for a DH potential) as a function of the density are 
shown again in figure ^ along with the run of the ratios 
Rd/Rb{Z = 4) and {rs)/RB{Z = 4). The Debye radius 
and the mean interparticle distance are calculated assuming 
X — 0.34, Y = 0.68 and Z = 0.02, a composition which is 
close to the one at the solar core today. 

We notice that when the density approaches the den- 
sity in the solar core, namely about 150 g cm~^, (a) the 
Debye radius becomes of the order of the mean interparticle 
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distance and hence the approximation of a smooth Debye 
screened potential loses its validity, emphasizing once more 
the conclusion reached in section ti. (b) In the solar core, we 
find that Rn ^ Rb{Z — 4) and therefore the probability for 
complete ionization of the Be is very high. 

However, the more important question here is the ratio 
of the mean interparticle distance to the Bohr radius since 
we are interested in the possibility that the ions of Beryllium 
still have bound electrons. Fig. m which depicts a graph for 
the value of (r^), indicates that this value at the center of 
the Sun it is close to Rb{Z = 4). Therefore, it is a delicate 
question whether the Beryllium ions possess any bound elec- 
trons. Finally, we point out that (rs), which is depicted in 
the figure, is the mean interparticle distance irrespective of 
their type. As we shall show, it is an underestimate in the 
case of a Beryllium ion embedded in Hydrogen and Helium 



4 PRESSURE IONIZATION AT T = 



take into account the close ion. As the speed of the perturb- 
ing ions is of the same order as the speed of the ion under 
consideration, the effect of the ions inside the Debye sphere 
cannot be averaged into a mean potential. The basic require- 
ment to consider the effect of the nearest neighbor directly 
(and not via a smooth potential averaged over many ions) 
leads to the idea of a Wigner-Seitz unit cell or ion-sphere. 

We can look also on the problem in the following way: 
If the bound state electronic wave function of a given ion 
overlaps significantly the bound electronic wave function of 
nearby ions the electron cannot be considered as bound (cf. 
Murillo fc Weisheit|[l998| ). 

We approached the problem of pressure ionization at 
T = in two steps which represent successive approxima- 
tions. In the first step we solve the Schrodinger equation 
for the Be'^ ion under the assumption that there is another 
nucleus at a distance (r^) away. In our particular case, the 
plasma contains ions with different charges and one cannot 
state that the ion sphere of all ions is identical. Hence the 
boundary condition must be imposed at (r^) = a{rs), where 
a is soon to be determined. In the pure periodic case one 



4.1 Tflit! bouiidtiry oontllllons 



shoul d app ly the Bloch condition ( Mardcr |2000 , see also 



et al 



Lai 



1991). We assume for simplicity spherical symmetry 



As discussed above, the classical method to evaluate the 
degree of ionization in a stellar plasma with a finite tem- 
perature is first to assume given mean distances between 
the particles, and assume that they are at rest, namely that 
T — 0. (However, we do keep the finite temperature in the 
calculation of the Debye radius). Once the energy levels are 
known, the effect of the temperature via the Boltzmann re- 
lation (leading to the Saha equation) is taken into account. 
In an actual plasma, the distance between the particles has 
a distribution and hence there is a distribution of cases. One 
assumes that the average of the results for the distribution 
is equal to the result for the average. We turn now to the 
T = case. 



and hence the Bloch condition becomes the requirement that 



(hp 
dr 



Th e que stion of pressure ionization is discussed by Chiu 



= 0. (11) 

r— a (rs ) 

The coefficient a is determined from the condition that the 
force vanishes at this point. Consequently, all wave func- 
tions we experimented with contained the condition that 
?/; vanishes at {rz) = ceifs), namely contained the factor 
(r — a{rs))^. Here {ts) is the distance to the nearest ion ir- 
respective of its charge (in spherical packing). Note that in 
the Sun Be^ is a trace element and hence the nearest neigh- 
bor would most probably be a proton or a Helium nucleus. 
The approximation can be considered as a muffin-tin poten- 
tial with a modified cell size, namely C oulomb inside (rz) 



&i Ng| ( |l999| ) within a general discuss ion about the energy '^^^ constant outside (cf. 



Lai et al 



1991). It is obvious that 



levels of ato ms in plasma a nd fo llows 



Roussel & O'Connell 



( [1974D and [Rogers et alj (|1970|), where the effect of the 
plasma was simulated by a screened potential. Pressure ion- 
ization depends primarily on the density. However, when the 
relevant scale is the Debye length, some effect of the tem- 
perature on the energy levels enters through the back door 
via the dependence of the Debye screening length on the 
temperature, which as stated before we do keep finite in the 
calculation of the potential (in the case that a DH potential 
in assumed). 

We distinguish between two possible situations: 

Case A: Rd > (rs) or No > 1 and 

Case B: Rd < (r,) or No < 1. 

The physical difference between the two cases is ex- 
pressed in the boundary conditions imposed on the wave 
function in the problem of the electronic structure of the 
ion. In the first case, the nearest neighbor is closer than the 
Debye distance and hence one expects it to affect the elec- 
tronic structure of the ion much more than the fact that the 
two ions are inside the same potential well. In the second 
case, the nearest neighbor is further away than the Debye 
radius and hence the boundary condition on the wave func- 
tion with a DH potential can be tl;{r oo) — 0. On the 
other hand, in the first case, the boundary condition should 



imposing the above condition on the trial function increases 
the eigenvalue and hence ionization would occur, assuming 
all other conditions are unchanged, at a lower density. 

Since we have a mixture of various ions, the point at 
which the force between ions vanishes varies with the ion. If 
our ion has a charge Z, then the corresponding ionic radius 
(rz) is given by: 



1 + 5 



(rs) where 5 — 



ZY^X.Z.jA, 



1/2 



(12) 



If all ions are equal, then {rz) = (l/2){rs). Clearly, the 
above definition is a generalization of the Wigner-Seitz cell 
idea to a mixture of species. In the present case, we assume 
the electron to be localized in the Wigner-Seitz cell. We 
will define complete ionization when the localization of the 
electron ceases. In the DH case, there is no such assumption. 
Nevertheless, there are other ions (and electrons) moving 
inside the Bohr radius. 

The relation of (rz) or (r^) to Rd is interesting. The 
condition (r^) = Rd can be written as 



pcrit 



(g cm-^) = 1.57(1 + 3X)^(^)' 



(13) 



For densities below the critical density, the mean interparti- 
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Figure 5. The run of the temperature (Tg), the density, and 
the critical density as well as the Debye radius and the mean 
interparticle distance (the latter given in units of the Bohr radius 
for a Z = 4 ion) in the Sun. 
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Figure 7. The energy of the ground state of an ions with charge 
Z immersed in a plasma inside a unit cell of size {rz)- 
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Figure 6. The temperature-density plane with the structure line 
of the Sun and the lines = {ts) for three different values of 
the Hydrogen mass fraction, 0.35, 0.5 and 0.7. The almost straight 
part of the structure line is the convective zone. 



cle distance is smaller than the Debye radius. Hence, when 
analyzing the possibility that an ion carries a bound elec- 
tron, the Debye length is not relevant. In fig. ^ we show the 
run of the density and temperature in the Sun along with the 
critical density (calculated for the actual temperature and 
composition). Fig. ^ depicts the domains of Rd = {ts) for 
three values of the hydrogen mass fraction, 0.35, 0.5 and 0.7. 
One finds that through most of the volume of the Sun, the 
density is always below the critical one and hence the anal- 
ysis of the structure of the electronic levels must take into 
account the nearby ion rather than the Debye radius. Only 
close to the surface does the situation change, the Debye ra- 
dius becomes smaller than the mean interparticle distance 
and the critical density becomes smaller than the actual den- 
sity. 



When one evaluates the pressure ionization for metals 
at r = 0, one assumes the Wigner Seitz cell. The rational 
for using it at higher temperatures is the fact that the speed 
of the electron in the bound state is so much greater than 
the speed of the ions, such that the ions can be assumed to 
be at rest. The us e of the ion-sphe re for opacity calculations 
was examined by Rozsnyai (1992). 



4.2 Schrodinger equation with Coulomb potential 

We used the variational principle method and trial func- 
tions of two typ es. The first type is taken from Roussel & 
O'Connel] (1974), namely a polynomial in r times an expo- 
nential function (the simple bound s-state), while the second 
type is a Fade approximation. 

The results for Hydrogen obtained using the two types 
of trial functions are compared and found to be practically 
the same to within a relative accuracy of 10~^ or better. 
Additional trials with other functions and parameters did 
not improve the results beyond the second significant digit. 

Interestingly, with our definition of {rz), the energy 
level is only a function of {rz)/RB[Z) and it is shown in 
fig. ^ Complete pressure ionization is found to occur at 
(rz) = 1.945i?B(Z) irrespective of Z. To obtain the results 
for a particular ion, one has to find its {rz) for the com- 
position, temperature and density under consideration. The 
results for T = are shown in fig. ^ We stress that these 
results are obtained for T = and do not depend on the 
Debye radius (and hence do not depend on the temperature 
indirectly). We conclude that neither Ne'^", nor species with 
a higher Z, are fully ionized in the solar core. 

The above results are easily translated into the condi- 
tions in the Sun. In fig. ^, we plot the run of the ground state 
binding energy of Be^ throughout the Sun. This calculation 
indicates that Beryllium is fully ionized in the Sun below a 
solar mass fraction of 0.66. 

The critical densities for the disappearance of the bound 
state in the corresponding Hydrogen like ions are given in 
Table ^. At a finite temperature, complete ionization takes 
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Figure 8. The translation of the previous figure to the density 
dependence of the energy levels of various species. All calculations 
assume that the Z 1 ions are trace elements. 
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Figure 9. The run of the ground level of Be^ in the Sun along 
with relevant thermodynamic parameters. On the left is the 
ground energy level axis. The axis for all other quantities is on 
the right. Temperature is given in 10*K, density in g cm~^. The 
arrow marks the point of complete ionization. 



places at somewhat lower densities because the temperature 
increases the excitation and with it the probability of ion- 
ization. 

We conclude that the densities at which full ionization 
of Be^ takes place are significantly lower than those found 
in the solar core. Thus, the small inaccuracies in the present 
estimate have no practical effect on the state of Be in the 
solar core. The effect on the entropy density of the envelope 
is yet another issue, and it will be discussed elsewhere. 

4.3 A Screened Potential? 

The Schrodinger equation was solved using the simple 
Coulomb potential. It did not incorporate the Debye 
screened potential. However, electron screening does take 



place and cannot be ignored. Should the Debye screened 
potential be used? Suppose that Rd ^ (ts) (which is not 
the case in the Sun), namely there are many particles in the 
Debye sphere. The Debye sphere contains electrons and ions, 
hence under these conditions there are ions which are closer 
to the specific ion than the Debye radius. These ions disturb 
the given ion and one should look for the bound level under 
these conditions, namely that there is another ion close by. 
This is exactly what has been done above. Thus, when the 
Debye radius is very large relative to the mean interparticle 
distance it creates a constant potential at the location where 
the two close particles are and the effect is a shift in the en- 
ergy and pressure ionization of the very high energy levels. 
On the other hand, when Rd ~ (rs) the Debye length and 
potential lose their meaning. This is exactly the situation in 
the solar core. A full treatment of this limit which incorpo- 
rate a Debye (or a more accurate) potential would have led 
to a still lower densities for the disappearance of the bound 
state). 



4.4 The Kohn-Sham equation 

There are two major deficiencies in the above treatment. 
The first one was discussed at length and has to do with the 
doubtful validity of the Debye potential under solar con- 
ditions and for the question of Be pressure ionization. A 
possible way t o overcome this problem is to us e a density 



functional (cf. Dharma-Wardana & Perrot 



1982 where it is 



applied to Hydrogen plasma and where bound states are 
found). The second deficiency is the neglect of the free elec- 
trons and their effect on the screening of the nucleus, correla- 
tions and exchange. We therefo re resorted to the Kohn-Sham 

in which these deficiencies 



equation (Kohn & Sham 



196; 



are taken care of, for finding the density at which pressure 
ionization takes place. The Kohn-Sham equation takes the 
A'^ electron wave function and treats it as a collection of sin- 
gle particle eigenfunctions. The governing equation of the 
Kohn-Sham density functional method is then: 



Ze 



r / r — r' 



(14) 



V'i(j-) = Eiilniv). 



where n(r) is the electron density and xpi is the single elec- 
tron wave function. The Kohn-Sh am equation ( hereafter KS 
equation) has additional merits (Mardei 2000). The major 
problem with Thomas Fermi equation is that high densities 
do not necessarily lead to high kinetic energies for the elec- 
trons. This problem is cured in the above KS equation. There 
are more merits of using the KS equation in our particu- 
lar problem. In the Hartree-Fock approach, the many-body 
wave function in form of a Slater determinant plays the key 
role in the theory. The Hartree-Fock equation if derived by 
minimization of the total energy is expressed by a determi- 
nant of wave functions which is extremely difficult to handle. 
In the density functional theory the key role is played by 
the observed quantity, the electron density. The Hohenberg- 
Kohn theorem then shows that for ground states the den- 
sity functional theory possess an exact energy functional and 
there exists a variational principle for the electron density. 
The KS equation is then an effective one electron equation 
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where the exchange operator in the Hartree-Fock equation is 
replaced by an exchange-correlation operator that depends 
only on the electron density. This is exactly what is needed 
in the present problem. The KS equation then treats the A'' 
electron problem as single electron wave functions. 

Let n(r) = X^iLi where the summation is car- 

ried over all electrons, be the electron density. In our par- 
ticular case, we examine the bound state of Be (Z=4) (as 
well as that of the higher Z trace species like C, N etc). The 
Be is a trace element immersed in a plasma of fully ionized 
Hydrogen and Helium (mostly) and negligible amounts of 
heavy elements. Hence, the major contribution to the elec- 
tron density comes from the electrons contributed by Hy- 
drogen and Helium. This term is essentially given by the 
environment in which the trace specie is immersed. Return- 
ing to the KS equation, the third term in the KS equation is 
the mutual electron-electron interaction between the bound 
electron and the free electrons which exist inside the 'effec- 
tive orbit' and it provides the effective electron screening of 
the ion. It is easy to estimate when this term becomes im- 
portant. The number of free electrons per Bohr radius, A''^ 
is given by; 



6.16 X IQ-^'^ 
Nb ^ ^3 



(15) 



This term becomes important for Nb ~ 1 or > 1.623 x 
IQ^i^s^ The forth term is the exchange which is given by 
dEex/dn where Eex is the exchange energy. Because of its 
unique properties, the KS equation gained popularity with 
physical chemists. The accuracy of using the KS equation 
for the calcu lation of ionization potentials in molecules is 
described in Curtiss et al. (1998). As a rule, the results of 
the KS equation are more accurate than those obtained from 
the Hartree-Fock approximation and reach the accuracy re- 
quired by quantum chemists. The implementation of the KS 
equation in a strophysics of d ense matter is described for the 
first time by 



Lai ct al 



(1991) 



We solved for the eigenvalue of the KS equation as- 
suming the composition of the present solar core. We used 
a variational principle with several trial functions since we 
are mostly interested in the eigenvalues and not in the wave 
functions. The results for Be^ are shown in figure |l^ along 
with the results for the Coulomb potential (and the same 
boundary condition). It appears that the correlations and 
exchange terms in the KS equation contribute to the fur- 
ther suppression of the energy level and complete ionization 
occurs a la KS at a higher density. Note that for sufficiently 
low densities, the KS predicts significant lowering of the 
ground state of a trace element relative to the continuum. 
This is a consequence of the exchange term which is mostly 
contributed by Hydrogen and Helium and not by the elec- 
trons of the trace ion under consideration. The phenomenon 
does not occur for species which are not trace elements (see 
later). The results with different trial functions vary a bit 
and shown are the best (lowest) results. The critical den- 
sities and (rz) as obtained in the two approximations are 
compared in Table |l|. From the table we see that in the 
Schrodinger approximation (rz) is constant while in the KS 
equation it varies and increases with the charge of the ion. 
A good fit for the value of the critical (rz) over this range 
of charges is rz/Z — 0.451 — 0.026^. This expression can be 
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Figure 10. A comparison between the MufBn-Tin models em- 
ploying the Schrodinger and the Kohn-Sham equations. Also 
shown is the radius of the atomic cell rz in units of the Bohr 
radius for a Hydrogen like ion with Z = 4 treated as a trace el- 
ement. The ground state energy level varies almost linearly with 
the density when close to the complete ionization density. 



Table 1. The Critical Density and the radius of the atomic cell 
for the Vanishing of a Bound State in Hydrogen Like Ions 



Ion 
Units 


Pcrit{Coul) 

g cm~^ 


(rz) 
Rb(Z) 


Pcrit{KS) 

g cm"'' 


(rz) 
Rb{Z) 


Be'' 


16.33 


1.921 


36.18 


1.474 


(jl2 


65.25 


1.923 


86.25 


1.752 


N14 


110 


1.922 


132.5 


1.806 


Ql6 


171 


1.927 


196.5 


1.840 


Mg24 


660 


1.925 


698.8 


1.890 



easily translated into a term added to the free energy so as 
to secure pressure ionization. 

Based on the KS equation we find that the CNO ele- 
ments are fully ionized in the core of the sun. Indeed, at 
r = the Oxygen still has a bound electron at the densities 
of the solar core, but the low binding energy and the high 
temperature impose complete ionization under the condi- 
tions in the center of the Sun. On the other hand, species 
with Z > 10, like Neon and Iron, still keep their K-shell 
electrons. Assuming that all species heavier or equal to Ne 
are fully ionized introduces a relative error of the order of 
<Z{1 - X{C) - X{N) - A"(O))/10, or about 10^=* in the 
pressure and speed of sound. 



4.5 The effect of the boundary condition in the 
Kohn-Sham equation 

The effect of the boundary condition on the result can be 
seen in the following way. We solved the KS equation under 
the condition il}{r — > oo) = for various densities. The re- 
sult is that as the density increases monotonically, the first 
bound state becomes monotonically more bound and pres- 
sure ionization never occurs. The free electrons are able to 
prevent ionization exactly like the Saha equation with elec- 
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tron degeneracy, which yields that the free electrons recom- 
bine as the density increases. 



4.6 Metallic Beryllium 



Ferrot ( |l990| ) applied the Neutral-Pseudo Atom (NPA) 
method to evaluate the equation of state and the degree 
of ionization of pure Be metal under normal conditions 
(p — 1.85 g cm~^ and T — Q) and at high densities, keeping 
T = 0. The results are not directly applicable to the Sun 
because Perrot discusses pure metallic Beryllium while we 
are interested in a trace Be atom embedded in a sea of pro- 
tons and a particles. Yet, the results are instructive for the 
comparison of the present method with others, especially 
because they serve as a consistenc y check. The applicability 
of the NPA method (Perrot 199C ) is limited t o compression 
ratios c below 40 (see fig. 5 in the Perrot 1990 and the expla- 



nation therein) while Be becomes fully ionized at compres- 
sion ratio of c = 50. At high compressions, band calculations 
show that the gap between the Is ba nd and the upper one 
close ( [Meyer- Ter-Vehn fc Zittel] |l98^ ) . The normal density 
of Beryllium is 1.85 g cm~^ and hence a compression ratio 
of 50 corresponds to a density of 92.5 g cm~^. Consequently, 
Metallic Beryllium is fu lly pres s ure io nized at T = and a 
density of 92.5 g cm" 



Perrot 



(199C) presents an extrapo- 
lation between the end point of his results (c = 40) and the 
band calculations (c = 50) . 

We note that when the density or temperature (or both) 
are increased starting from p — 92.5 g cm~^ and T = no 
bound state can re-appear. There is simply no bound state 
at higher densities and/or higher temperatures if it does not 
exist for T = and a given density. 

Our calculations of the pressure ionization assume that 
the heavy elements under considerations are trace elements. 
Yet, it is of interest to compare the present method with oth- 
ers. The above results for metallic Beryllium provide such an 
opportunity. When the specie under consideration is a trace 
element, the electrons are contributed by the Hydrogen and 
Helium and there is no connection between the number of 
electrons in a unit cell and the charge of the ion. (There is 
no condition of charge neutrality per each ion sphere. The 
number of electrons in the Beryllium ion sphere is not neces- 
sarily 4). In the pure Beryllium case, when we search for the 
density of complete ionization, we assume that the free elec- 
trons are the the first three electrons of Beryllium and that 
the forth one is bound. Thus the exchange term is evaluated 
on the basis of Beryllium ionized three times. 

As can be deduced from the previous calculation, a crit- 
ical factor is the packing of the specie. If spherical packing 
is assumed, then the KS equation predicts co mplete ioniza- 

(|199C|) 



Perrot 



tion at p ~ 56 g cm~ . On the other hand 
quotes that complete ionization is reached at 92.5 g cm~^. 
The exact lattice structure of Beryllium at very high densi- 
ties is not known and Be apparently undergoes a structural 
change at a compression ratio of about 3. However, if we 
as sume an f e e latt ice structure and resort to the appendix 
in Lai et al, (1991) to find the relation between the lattice 
structure and the radius of the atomic cell (r^ = (\/2/4)a) 
we find complete ioniz ation a t a den sity of 89.3 g cm~^, in 
good agreement with Perrot (199C). This result reassures 



5 THE EFFECT OF BE FULL IONIZATION ON 
THE SOLAR NEUTRINO FLUX 

When Beryllium is completely pressure ionized in the core 
of the Sun, electron capture by Be takes place via the con- 
tinuum only. Consequently, no corrections to the rate due 
to bound electrons should apply. 

The effect of the complete ionization of Be^ on the solar 
neutrino flux can be easily estimated without recourse to 
detailed solar models because Be'^ is a trace element and the 
amount of energy released by its reactions is negligible. Two 
facts are relevant, (a) The Be^ electron capture neutrino flux 
and the B* flux are proportional to the abundance of Be^. 
(b) The rate of electron capture is much larger than the 
rate of proton capture. Consequently, a decrease of about 
20-30% in the rate of the Be^ -|- e~ capture increases the 
abundance of Be^ by the same factor and hence increases the 
B* neutrino flux by the same factor. The flux of the Be'^-l-e" 
neutrinos is unchanged because the decrease in the rate is 
fully compensated for by the increase in the abundance. As a 



matter of fact, the rate of Be +e 
condition that the rate of Be^ -I- e 



is determined (under the 
is much larger than the 
rate of Be^ + p ) only by the rate of He^ -|- He*. The above 
treatment of the pressure ionization of Be'' applies to all 
solar models irrespective of any other parameter. The effect 
on the predicted fiux assuming few different sets of nuclear 
parameters is shown in Table |^ The four models referred 



to in t he ta ble are the follow ing. Cas97 ref e rs to Castellani 
(|1997|), DSm refers to [Par fc Shavi^ (|l996[), RMP98 



et al 



refers to 



Adelberger fc et al. 



(1998) and NACRE refers to 



Angulo & et al. (1999[). The models differ mainly in the 



that our analysis for Be*^ as a trace element is valid as well 



parameters for the rates of the nuclear reaction and slightly 
in the details of the equation of state. 



6 CONCLUSION 

In this work, we were interested in the internal structure of 
a trace ion in a high density plasma, and not in the gen- 
eral properties of the plasma. The relevant physical quan- 
tity for the latter, is the effective charge of an ion — not the 
electronic structure provided the effective radius of the elec- 
tronic wave function is much smaller than the Debye radius. 
As to the structure of the ion itself, we distinguished between 
two cases depending on whether the mean interparticle dis- 
tance is smaller or larger than the Debye radius. We found 
that throughout most of the Sun, the mean interparticle dis- 
tance is slightly smaller than the Debye radius. 

Under such conditions, the ionization state of a given 
ion depends mostly on the distance to the nearest neigh- 
bor rather than the distance scale of the screened potential. 
Hence, we approximate the condition of the plasma by a gen- 
eralized Wigner-Seitz cell around each ion plus the proper 
boundary condition on the surface of the cell. We then im- 
pose the Bloch condition on the boundary of the generalized 
Wigner-Seitz cell. The fact that the nearest neighbor is so 
close implies that the implementation of this condition has a 
major effect on the point at which complete ionization takes 
place. 

We have calculated the critical density of complete ion- 
ization assuming at first the Schrodinger equation and a 
Coulomb potential. Next, we implemented the Kohn-Sham 
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Table 2. The effect of treating the Be as fully ionized on the solar neutrino flux 



Cas97 DS96 RMP98 NACRE 





full 


partial 


full 


partial 


full 


partial 


full 


partial 


0(S8)/io6cm-2s-l 


6.36 


5.33 


3.74 


3.13 


6.20 


5.19 


6.30 


5.28 


0(S8)loV<^(Se'^) 


1.42 


1.19 


0.935 


0.783 


1.25 


1.05 


1.25 


1.04 


Ga(SNU) 


133 


133 


121 


121 


132 


132 


137 


137 


C1{SNU) 


9.14 


7.80 


5.90 


5.10 


8.96 


7.65 


9.29 


7.96 



equation which takes into account the exchange interaction 
due to the free electrons and the screening by them. Com- 
parison of our results with the NPA method yields a very 
good agreement for the case of metallic Beryllium. 

We mention that when the boundary condition imple- 
mented in the KS equation is i/;(r — > oo) — the phe- 
nomenon of pressure ionization disappears. 

The improvements in the treatment of pressure ioniza- 
tion of trace elements show that all species with Z < 8 are 
fully ionized in the core of the Sun. (The critical density for 
Oxygen is a bit higher than the density in the core of the Sun 
but the high temperature secure the complete ionization of 
Oxygen). The frequently applied correction to the Be^ elec- 
tron capture rate due to bound electrons, is not needed. The 
revised prediction for the (^^(B*) from the Sun is higher by 
about 20-30%. 

The discussion here was centered around the high den- 
sity effects. Howev er, one should take into ac count collective 
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effects as well (cf. Murillo & Weisheit 
the electron plasma oscillations is 



1998). The energy of 



hcje = 3.7 X 10' 



eV. 



(16) 



At an electron number density of 10 the energy of the 
plasma oscillations amounts to 370 eV more than the ion- 
ization potential of Be'^. Thus the bound states of Be^ have 
been broadened into continuum states long before. 
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